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1.  Introduction 


Various  2-azidoethanamines  (EAs)  have  demonstrated  potential  as  fuels  for  hypergolic  rocket 
propulsion  systems  (Thompson  et  ah,  1998;  Stevenson  et  ah,  2004;  McQuaid,  2006),  and  as  part 
of  an  effort  to  determine  whether  candidates  from  this  class  warrant  further  testing  and 
development,  data  that  can  be  used  to  assess  the  risks  they  pose  to  human  health  and  the 
environment  are  being  generated  and  reported  with  support  from  an  Environmental  Quality 
Technology  (EQT)  program.  The  data  that  has  been  assembled  to  date  suggest  that  the  class  (in 
general)  is  not  as  hazardous  as  the  hydrazine-based  fuels  that  are  currently  standard  in  hypergolic 
propulsion  systems  (Bausum  and  Johnson,  2006).  However,  much  of  the  data  for  individual 
candidates  has  been  estimated,  and  due  to  a  paucity  of  relevant  measured  values  in  the  open 
literature  (and  the  EQT  program’s  limited  ability  to  obtain  them  otherwise),  some  estimates  are 
based  on  methods  that  are  not  well  validated.  As  discussed  in  a  report  summarizing  a  prior  EQT- 
sponsored  study  to  establish  physical  property  data  for  the  class  (McQuaid,  2008),  estimates  for 
Henry’s  law  constants  are  among  this  group. 

Representing  an  expectation  for  the  ratio  of  a  compound’s  concentrations  that  will  exist  in  the 
two  phases  of  a  ternary  compound-air-water  mixture  at  equilibrium,  a  compound’s  Henry’s  law 
constant  is  employed  to  predict  its  partitioning  between  air  and  water  in  the  environment. 
Expressible  in  various  ways,  the  values  of  the  constants  that  are  presented  in  this  report  are 
defined  per 

HLC(7)  =  RT\cair(T)lcwateJ{T)\,  ( 1) 

where  R  is  the  universal  gas  constant  (8.23  x  10  atm-m  /mole-K),  T  is  the  temperature  (in  K), 
cai,-(T)  is  the  concentration  of  the  compound  in  the  mixture’s  gaseous  phase,  and  cwate>(  T)  is  the 
concentration  of  the  compound  in  the  mixture’s  aqueous  phase.  In  addition,  for  environmental 
risk  assessment,  it  is  typical  to  seek  and/or  establish  HLC(7)  values  for  mixtures  at  298  K  and 
1  atmosphere  (atm),  and  HLC(7)  values  at  that  temperature  and  pressure  are  the  only  ones 
considered  here. 

In  the  aforementioned  EQT-sponsored  study,  an  attempt  was  made  to  establish  the  reliability  of  a 
quantitative  structure  property  relationship’s  (QSPR’s)  HLC(298)  estimates  for  EAs. 
Implemented  by  the  Environmental  Protection  Agency’s  (EPA’s)  HENRYWIN  program 
(HENRYWIN,  2003),  the  QSPR  is  a  bond  contribution  model,  and  it  has  a  coefficient  that  was 
specifically  developed  for  azide  groups.  It  was  surmised  that  the  coefficient  was  derived  from  a 
measured  HLC(298)  value  for  hydrazoic  acid  (HN3)  (McQuaid,  2008),  but  the  only  check  that 
was  able  to  be  made  of  the  QSPR’s  validity  for  EAs  was  a  comparison  of  its  estimate  for 
2-azido-N,N-dimethylethanamine  (DMAZ)  to  a  value  derived  from  “separately  measured” 
values  for  DMAZ’s  vapor  pressure  [VP(I)]  and  water  solubility  [WS(7)]  at  298  K, 
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HLC(298)  ~  VP(298)/WS(298). 


(2) 


(That  is,  the  value  was  derived  from  [1]  a  VP(298)  value  measured  above  neat  liquid  DMAZ 
[Dee,  2000]  and  [2]  a  separate  WS(298)  measurement  [Manzara,  2007],  not  the  measured 
concentrations  of  DMAZ  in  the  two  phases  of  a  ternary  DMAZ-air- water  mixture  at 
equilibrium.)  The  value  found  via  the  equation  2  approximation  (9.7  x  10  atm-m  /mol)  is 
almost  2  orders  of  magnitude  higher  than  the  HENRYWIN  estimate  (1.2  x  10  atm-m  /mol). 

Since  the  identification  of  an  HLC(298)  estimation  method  whose  predictions  for  EAs  could  be 
trusted  to  be  within  an  order  of  magnitude  of  actual  values  was  desired,  the  aforementioned 
comparison  was  not  encouraging.  However,  two  other  considerations  suggested  that  the  QSPR’s 
estimate  for  DMAZ  might  be  reasonable.  The  first  was  that  an  HLC(298)  value  for  HN3  derived 
from  equation  2  with  a  measured  WS(298)  value  and  a  trusted  VP(298)  estimate  was  also  2 
orders  of  magnitude  higher  than  a  “true”  HLC(298)  measurement  for  the  compound.  Thus,  it 
could  be  imagined  that  the  HLC(298)  estimate  derived  for  DMAZ  from  equation  2  was  too  high, 
a  speculation  being  that  hydrogen  bonding  between  gas-phase  water  molecules  and  DMAZ’s 
lone  pair  sites  leads  to  a  lower  DMAZ  concentration  above  water  than  one  would  expect  on  the 
basis  of  its  concentration  above  its  neat  liquid  form.  (A  DMAZ  molecule  has  hydrogen-bond 
accepting  lone  pair  sites  that  allow  it  to  hydrogen  bond  with  water,  but  it  does  not  have  a 
hydrogen-bond  donating  site  that  allows  it  to  hydrogen  bond  to  another  DMAZ  molecule.) 

Secondly,  HENRYWIN’ s  HLC(298)  estimate  for  DMAZ  could  be  rationalized  on  the  basis  of  its 
comparison  to  a  relatively  well-validated  estimate  for  N,N,N’,N’-tetramethylethylenediamine 
[TMEDA,  (CH3)2NCH2CH2N(CH3)2,  2.4  x  1CT8  atm-m3/mol].  Both  DMAZ  and  TMEDA  can  be 
thought  of  as  N,N-dimethylethanamine  with  a  functional  group  at  the  2-position  of  the  ethyl 
chain;  in  DMAZ’s  case  the  functional  group  being  an  azide,  and  in  TMEDA’ s  case  the  functional 
group  being  a  dimethylamino  fragment.  Therefore,  one  can  predict  with  some  confidence  the 
relative  differences  in  their  VP(298)  and  WS(298)  values  and,  by  extension,  their  HLC(298) 
values.  TMEDA  and  DMAZ  have  similar  VP(293)  values;  i.e.,  10  mm  Hg  (Acros  Organics, 
1996)  vs.  8  mm  Hg  (Dee,  2000),  which  is  an  observation  one  would  predict  on  the  basis  of  their 
molecules  having  similar  molecular  weights  and  little  potential  for  strong  attractive  interactions 
(such  as  hydrogen  bonding)  with  another  of  their  own  kind.  But  TMEDA’ s  (second) 
dimethylamino  fragment  can  bond  more  strongly  to  water  than  DMAZ’s  azide  group  can,  which 
leads  to  the  expectation  that  TMEDA  will  have  a  higher  WS(298)  than  DMAZ.  A  value  for 
DMAZ’s  WS(298)  that  was  derived  from  measurements  reported  by  Manzara  (2007)  indicates 
that  it  is  ~1.7  x  105  mg/L  while  TMEDA’s  WS(298)  is  estimated  to  be  “1  x  106  mg/L” 

(McQuaid,  2008).  (A  compound  with  a  WS(298)  value  of  1  x  10(1  mg/L  in  essence  means  that  it 
is  miscible  with  water  in  all  proportions.)  Thus,  while  quantitative  predictions  based  on  these 
observations  are  not  strongly  defensible,  the  prediction  that  DMAZ’s  HLC(298)  would  be  a 
factor  of  5  higher  value  than  TMEDA’s  seemed  reasonable.  Thus,  HENRYWIN ’s  estimates  for 
EA  candidates  were  considered  adequate  for  the  purpose  of  risk  assessment  for  the  stage  of  the 
selection  process  that  existed  at  the  time,  and  they  were  reported  (McQuaid,  2008). 
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This  report  considers  HLC(298)  estimates  for  several  EA  candidates  that  were  obtained  by  an 
alternate  method.  Investigated  after  the  completion  of  the  prior  study,  the  method  is  based  on 
calculating  standard-state  (298  K  and  1  atm)  gas-  and  aqueous-phase  Gibbs  free  energies: 

G°g  (298)  and  G°  (298)  ,  respectively,  from  results  obtained  with  a  density  functional 

theory/polarizable  continuum  model  (DFT/PCM).  From  G°g  (298)  and  G°q(298) ,  one  can 
compute  a  compound’s  standard-state  Gibbs  free  energy  of  solvation  [  AG°0,(298)  ], 

AG.:/(298)  =  g;(298)-G;(298),  (3) 

and  it  is  related  to  a  compound’s  HLC(298)  per 

HLC(298)  =  (R  *  298)  exp[AG“,  (298)  !{R  *  298)] .  (4) 

Having  more  bona  fide  theoretical  foundations  than  HENRYWIN’s  QSPR,  DFT/PCM-based 
methods  can  in  certain  cases  be  a  preferable  option  for  generating  HLC(298)  estimates. 

However,  for  compounds  whose  molecules  have  a  relatively  large  number  of  atoms,  and  thus  a 
relatively  large  number  of  molecular  conformations  that  need  to  be  characterized  in  order  to 
compute  AG", (298)  ,  the  method  is  far  more  expensive  to  implement  than  a  QSPR.  Moreover, 

such  methods  have  their  own  potential  pitfalls,  and  their  results  should  be  validated.  (Each 
1-kcal/mol  error  in  a  AGS°0,(298)  calculation  will  produce  over  a  factor-of-5  error  in  the 

HLC(298)  estimate,  and  errors  of  such  magnitude  are  about  the  smallest  one  can  hope  to  achieve 
with  such  methods.)  But  the  lack  of  relevant  measured  values  needed  to  validate  EA  HLC(298) 
estimates  has  already  been  noted.  Thus,  no  attempt  was  made  to  employ  such  a  method  for  the 
prior  study.  However,  DFT/PCM  results  generated  for  a  recently  completed  study  of  the 
aqueous-phase  dissociation  constants  ( Ka )  of  the  conjugate  acids  of  DMAZ,  2-azido-N- 
methylathanamine  (MMAZ),  l-(2-azidoethyl)pyrrolidine  (PyrrAZ),  and  4-(2-azidoethyl)- 
morpholine  (MorphAZ)  (Stevenson  and  McQuaid,  2009)  presented  a  low-risk,  low-cost 
opportunity  to  obtain  DFT/PCM-based  HLC(298)  estimates  for  them,  and  it  was  taken. 

This  report  describes  the  DFT/PCM-based  method  employed  to  derive  HLC(298)  estimates  and 
presents  the  method’s  estimates  for  DMAZ,  MMAZ,  PyrrAZ,  MorphAZ,  dimethylamine 
(DMA),  trimethylamine  (TMA),  and  1-methylpyrrolidine  (MP).  (The  chemical  structures  of  the 
EAs  are  depicted  in  figure  1 .)  Inferences  that  can  be  drawn  from  a  comparison  of  the  estimates 
to  values  obtained  by  other  means  are  discussed. 
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MMAZ  DMAZ  PyrrAZ  MorphAZ 


Figure  1.  The  chemical  structures  of  EAs  discussed  in  this  report. 


2.  Computational  Method 


The  theoretical  method  employed  to  calculate  G°(298)  and  G"  (298)  values — and  from  them 

AG"o/(298)  and  HLC(298) — is  based  on  calculating  gas-  and  aqueous-phase  Gibbs  free  energies 

for  equilibrium  conformers  of  molecules  of  the  compound  of  interest.  Denoted  for  a  given 
conformer  (/)  as  G„  (298,/)  and  Gaq  (298,  i),  respectively,  figure  2  depicts  the  relationship 

between  these  two  parameters  and  their  reference  state,  which  for  a  given  molecule  is  the 
conformer  (hereafter  designated  by  the  symbol  X)  with  the  lowest  Gg  (298,  /)  value.  Gaq  (298,  /) 

and  Gg  (298,  X )  are  connected  by  the  relationship 

Gaq  (298,  /)  =  Gg  (298,  X)  +  AGg  (298,  /)  +  AEg  (/)  +  AGsol  (298,  /),  (5) 

where  (1)  AGa(298,/)  is  the  difference  between  the  free  energy  of  the  gas-phase  equilibrium 

conformer  whose  geometry  most  closely  resembles  the  aqueous-phase  conformer  and  the  free 
energy  of  the  lowest-energy  gas-phase  conformer  [  Gg  (298,  /)  -  G„  (298,  X)  ],  (2)  AG,  (/)  is  the 

difference  between  the  gas-phase  electronic  energy  of  the  conformer’ s  aqueous-phase  geometry 
and  the  electronic  energy  of  the  gas-phase  conformer  whose  geometry  it  most  closely  resembles, 
and  (3)  AGvo/  (298,/)  is  the  free  energy  change  of  solvation  for  the  aqueous-phase  geometry. 

Gg(298,/)  is  computed  from 

G„  (298,  /)  =  Eg  (/)  +  ZPVEg  (/)  +  A  Gcgorr  (298,  /),  (6) 

where  E  (i)  is  the  conformer’s  electronic  energy  at  0  K,  ZPVEg(i )  is  its  zero-point  vibrational 
energy,  and  AG™"'  (298,/)  is  its  “ PV-TS ”  energy  at  298  K,  where  P,  V,  and  S  are  pressure, 
volume,  and  entropy,  respectively. 
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Values  for  the  parameters  defined  in  the  preceding  paragraph  were  derived  from  a  slightly 
modified  version  of  a  procedure  outlined  by  Saracino  et  al.  (2003)  and  calculations  performed 
with  the  Gaussian  03  (G03)  suite  of  quantum  chemistry  codes  (Frisch  et  ah,  2003).  Geometry 
optimizations  were  performed  to  identify  the  equilibrium  structures  that  the  molecules  of  interest 
would  acquire  in  an  aqueous  environment.  A  PBEPBE/6-31+G(d,p)  model  (Perdew  et  al.,  1996, 
1997;  Petersson  et  al.,  1988;  Clark  et  al.,  1983)  coupled  with  the  PCM  developed  by  Cossi  et  al. 
(2002)  was  employed  for  these  optimizations.  The  optimizations’  convergence  criteria  were 
maximum  force  <0.000450  hartree/bohr,  root  mean  square  (RMS)  force  <0.000300  hartree/bohr, 
maximum  displacement  <0.001800  bohr,  and  RMS  displacement  <0.001200  bohr.  In 
accordance  with  the  procedure  outlined  by  Saracino  et  al.,  each  structure’s  A Gsoi(298,i)  value 
was  determined  from  the  difference  between  its  HF/6-31+G(d,p)-calculated  gas-  and  (PCM- 
simulated)  aqueous-phase  single  point  electronic  energies.  “UAFIF”  radii  were  specified  for  the 
PCM  calculations  (Barone  et  al.,  1997). 

The  geometries  of  the  equilibrium  structures  found  with  the  PCM-mediated  PBEPBE/6- 
31+G(d,p)  model  were  subsequently  employed  as  starting  points  for  gas-phase  PBEPBE/6- 
3 1+G(d,p)  geometry  optimizations.  The  same  convergence  criteria  as  those  employed  for  the 
PCM-mediated  optimizations  were  employed  for  these  optimizations.  For  all  structures  so 
obtained,  normal  mode  analyses  were  conducted,  and  ZPVEg(i )  and  AG™"' (298,/)  values  were 

derived  from  them.  G03’s  calculations  of  AG™"  (298,/)  are  based  on  the  ideal  gas  law  and  the 
calculation  of  S  via  a  standard  method  (McQuarrie,  1999). 
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An  issue  in  the  calculation  of  AG™"'  (298,/)  that  was  considered  is  the  fact  that  many  of  the 

molecules  characterized  have  low-frequency  vibrational  modes  whose  energy  levels  (and 
therefore  -TS  contributions  to  the  free  energy)  are  more  properly  described  by  hindered  internal 
rotation  formulae  than  by  harmonic  oscillator  formulae.  Therefore,  the  G03  option  for  applying 
the  more  appropriate  formulae  was  attempted  for  all  cases.  However,  the  method  failed  for  all 
the  ring-containing  structures.  The  default  G03  method,  which  is  to  treat  all  vibrational  modes 
as  harmonic  oscillators,  was  therefore  employed  for  these  cases.  It  was  assumed  that  the 
difference  in  the  results  produced  by  the  two  procedures  would  be  small,  and  that  expectation 
was  met  in  the  (only)  case  examined  to  test  the  assumption.  The  case  was  DMAZ,  and  for  all  the 
conformers  that  have  been  identified  for  it,  the  default  calculation  and  the  hindered  rotor  option 
yielded  AG™"  (298,/)  values  that  were  identical  to  10'6  hartrees.  Therefore,  this  matter  was  not 

further  considered. 


It  should  also  be  noted  that  some  of  the  gas-phase  geometry  optimizations  produced  structures 
that  were  not  geometrically  analogous  to  the  aqueous-phase  equilibrium  geometries  from  which 
the  optimizations  were  started.  Regardless  of  the  reason,  if  a  gas-phase  geometry  optimization 
bypassed  a  local  minimum  that  closely  resembles  the  aqueous-phase  geometry  (starting  point) 
and  ended  in  a  lower  energy  minimum  that  does  not,  the  sum  of  A Gg (298,/)  and  A Eg(i)  that  is 

obtained  is  still  the  same  as  it  would  be  had  the  “optimal”  local  minimum  been  found.  Only  the 
relative  contributions  of  the  two  parameters  to  the  total  are  different.  Thus,  from  the  standpoint 
of  having  the  potential  to  be  sources  of  error  in  the  calculation  of  the  observable,  such  cases  are 
not  a  concern. 

To  complete  the  calculation  of  parameters  employed  to  compute  Ga?(298,/)  values,  gas-phase 

PBEPBE/6-31 1+G(2d,2p)  single  point  energies  for  the  gas-  and  PCM-mediated  PBEPBE/6- 
3 l+G(d,p)-detennined  geometries  were  obtained.  For  each  conformer,  A Eg(i)  was  computed 

from  the  difference  between  the  two  results,  and  the  result  for  the  gas-phase  structure  was  added 
to  the  structure’s  (unsealed)  ZPVEg{i)  and  AG™"  (298,/)  values  to  obtain  Gg(298,z)  values,  and 

from  them  A Gg  (298,  /)  values.  The  Gaq  (298,  /)  values  were  obtained  by  adding  A Gg  (298,  /) , 
AEg  (/)  ,  and  AGso/(298,/) . 


To  compute  G°(298)  and  G°  { 298)  for  a  molecule’s  ensemble  of  conformers,  it  was  assumed 

that  the  relative  populations  of  the  conformers  in  each  phase  would  conform  to  a  Boltzmann 
distribution, 


Paqlg( 298,/)  = 


g(/)exp(-AGa?/g(298,/)/R*298) 


£g(/)exp(-AGat?/g(298,/)/R*298) 


(V) 
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where  g(i)  is  conformer/’s  degeneracy  factor.  G°  (298)  and  G°  ( 298)  were  then  calculated  by 
weighting  the  sums  of  G',,  (298,/)  and  Ga</(298,z)  in  accordance  with  those  distributions,  i.e., 

G;;s  (298)  =  2  Pa„g  (298,  i)Gatlg  (298 ,i).  (8) 

i 

In  the  implementation  of  equation  7,  all  confonners  of  all  molecules  were  assumed  to  have  C/ 
symmetry  and  therefore  assigned  a  g(i)  value  equal  to  1 . 


3.  Results  and  Discussion 


Table  1  lists  the  D F T/ P C M - b a s e d  A G °ol  (298)  and  HLC(298)  estimates  that  were  obtained  for 

DMA,  TMA,  MP,  DMAZ,  MMAZ,  PyrrAZ,  and  MorphAZ,  and  it  compares  the  FILC(298) 
estimates  to  values  obtained  by  other  means.  The  DFT/PCM  results  employed  as  input  for  the 
calculations  are  given  in  the  appendix. 

Table  1.  Measured  and  estimated  HLC(298)  values  for  various  amines. 


Compound 

A  °Gsnl 

HLC(298) 

HLC(298) 

HLC(298) 

DFT/PCM 

HENRYWIN 

Measured 

(kcal/mol) 

(atmm3/mol) 

(atmm3/mol) 

(atmm3/mol) 

DMAa 

-5.3 

3.8  x  10~6 

1.7  x  10~5 

1.8  x  10~5 

TMAa 

-3.5 

6.6  x  10~5 

3.6  x  10~5 

1  x  10^ 

MPb 

-4.5 

1.2  x  10~5 

2.8  x  10~5 

3  x  lO"5 

DMAZ 

-7.6 

6.4  x  10~8 

1.2  x  10-7 

— 

MMAZ 

-10.2 

2.9  x  10~10 

5.3  x  10~8 

— 

PyrrAZ 

-9.2 

4.3  x  10~9 

9.1  x  10~8 

— 

MorphAZ 

-13.2 

5.1  x  10~12 

8.0  x  10~10 

— 

“Christie  and  Crisp  (1967) 
bCabanietal.  (1971) 


Because  measured  HLC(298)  values  were  found  for  DMA,  TMA,  and  MP,  the  results  for  these 
cases  represent  a  test  of  the  reliability  of  the  DFT/PCM-based  method.  (The  measured  values 
were  in  the  training  set  from  which  FIENRYWIN’s  QSPR  was  derived,  so  they  do  not  constitute 
validating  data  for  it.)  The  DFT/PCM-based  HLC(298)  estimate  for  DMA  (3.8  x  10  6 
atm-rn  /mol)  is  a  little  less  than  a  factor  of  5  lower  than  the  identified  measured  value  (1.8  x  10 
atm-m3/mol),  the  DFT/PCM-based  HLC(298)  estimate  for  TMA  (6.6  x  10  5  atm-m3/mol)  is  a 
little  less  than  a  factor  of  2  lower  than  the  identified  measured  value  (1  x  10  4  atm-m3/mol),  and 
the  DFT/PCM-based  HLC(298)  estimate  for  MP  (1.2  x  10  atm-m3/mol)  is  a  little  over  a  factor 

_ c  o 

of  2  lower  than  the  identified  measured  value  (3  x  10  atm-m  /mol).  Combined,  the  results  for 
DMA,  TMA,  and  MP  suggest  that  the  DFT/PCM  method  is  valid  for  aliphatic  amines,  but  that 
its  predictions  for  compounds  in  the  class  will  be  somewhat  low. 
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The  DFT/PCM-based  HLC(298)  estimate  for  DMAZ  (6.4  x  10  8  atm-m3/mol)  is  about  a  factor  of 
2  smaller  than  the  HENRYWIN  estimate.  As  such,  the  result  (coupled  with  the  results  for  DMA, 
TMA,  and  MP)  is  considered  support  for  (1)  the  contention  that  HENRYWIN’s  HLC(298) 
estimate  for  DMAZ  is  reasonable  and  (2)  the  suspicion  that  the  HLC(298)  value  derived  for 
DMAZ  from  VP(298)  and  WS(298)  values  is  too  high. 

The  DFT/PCM-based  HLC(298)  estimate  for  MMAZ  (2.9  x  10  10  atm-m3/mol)  is  more  than  2 

_o  o 

orders  of  magnitude  less  than  the  HENRYWIN  estimate  (5.3  x  10  atm-m  /mol).  Similarly,  the 
DFT/PCM-based  HLC(298)  estimate  for  PyrrAZ  (4.3  x  10  9  atm-m3/mol)  is  over  an  order  of 
magnitude  less  than  the  HENRYWIN  estimate  (9.1  x  10  s  atm-m3/mol),  and  the  DFT/PCM- 
based  HLC(298)  estimate  for  MorphAZ  (5.1  x  10  '  atm-m  /mol)  is  over  2  orders  of  magnitude 
less  than  the  HENRYWIN  estimate  (8  x  10  10  atm-m3/mol).  Thus,  one  might  conclude  that  the 
agreement  between  the  two  methods’  estimates  for  these  three  compounds  is  poor.  However,  as 
has  been  noted,  the  volatilization  from  water  of  compounds  with  HLC(298)  values  below  3  x 

_n  •j 

1 0  atm-m  /mol  is  typically  assumed  to  be  negligible.  Therefore,  rather  than  being  significantly 
different,  the  estimates  from  both  the  DFT/PCM-based  method  and  HENRYWIN  lead  to  this 
conclusion  for  these  three  compounds. 


4.  Summary 


HLC(298)  estimates  for  four  EAs  and  three  aliphatic  amines  were  derived  from  a  DFT/PCM- 
based  method.  The  four  EAs  are  2-azido-N,N-dimethylethanamine  (DMAZ),  2-azido-N- 
methylethanamine  (MMAZ),  l-(2-azidoethyl)pyrrolidine  (PyrrAZ),  and  4-(2- 
azidoethyl)morpholine  (MorphAZ),  and  the  three  aliphatic  amines  are  dimethylamine  (DMA), 
trimethylamine  (TMA),  and  1-methylpyrrolidine  (MP).  The  DFT/PCM-based  estimates  for  the 
three  aliphatic  amines  indicate  that  the  method’s  estimates  for  aliphatic  amines  will  be  within  an 
order  of  magnitude  of  actual  values,  but  somewhat  low.  The  DFT/PCM-based  HLC(298) 
estimate  for  DMAZ  is  about  a  factor  of  2  smaller  than  the  estimate  obtained  for  it  from  the 
EPA’s  HENRYWIN  program  and  thus  supports  the  reasonableness  of  the  HENRYWIN  estimate. 
In  percentage  terms,  the  DFT/PCM-based  HLC(298)  estimates  for  MMAZ,  PyrrAZ,  and 
MorphAZ  are  all  significantly  lower  than  those  obtained  from  HENRYWIN,  but  the  estimates 
from  both  methods  suggest  that  the  volatilization  of  these  compounds  from  water  will  be 
negligible.  Taken  together,  the  results  from  the  DFT/PCM-based  method  support  the  conclusion 
of  a  previous  study  of  HENRYWIN’s  HLC(298)  predictions  for  EA  candidates  (McQuaid, 

2008),  i.e.,  that  the  predictions  merit  use  as  input  for  modeling  such  candidates’  transport 
through  the  environment. 
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Appendix.  DFT/PCM  Results  for  Various  Amines 


This  appendix  reproduces  a  subset  of  density  functional  theory/polarizable  continuum  model 
(DFT/PCM)  results  that  were  generated  for  a  study  that  was  performed  to  determine  if  the 
aqueous  phase  dissociation  constants  ( Ka )  of  the  conjugate  acids  of  2-azidoethanamines  (EAs) 
correlated  with  features  of  the  molecular  structure  of  the  acids  or  their  corresponding  bases. 
Included  here  are  all  the  results  that  are  needed  to  compute  AG", (298)  values  for 

dimethylamine  (DMA),  trimethylamine  (TMA),  1-methylpyrrolidine  (MP),  2-azido-N,N- 
dimethyl  ethanamine  (DMAZ),  2-azido-N-methylethanamine  (MMAZ),  l-(2-azidoethyl)- 
pyrrolidine  (PyrrAZ),  and  4-(2-azidoethyl)morpholine  (MorphAZ).  For  DMA,  TMA,  and  MP, 
only  one  conformer  is  important.  The  molecular  structures  of  DMAZ’s  and  MMAZ’s 
conformers  have  been  discussed  in  detail  in  prior  reports,1,2  and  their  labels  in  the  tables  in  this 
appendix  correspond  to  the  labels  in  those  reports.  The  labels  of  the  conformers  of  PyrrAZ  and 
MorphAZ  were  assigned  on  the  basis  of  their  molecular  structure,  a  given  letter  identifying  a 
correspondence  between  the  molecular  geometry  of  a  PyrrAZ  or  MorphAZ  conformer  and  that 
of  an  MMAZ  or  DMAZ  conformer  with  the  same  letter. 


Table  A-l.  DFT/PCM  results  for  DMA,  TMA,  and  MP. 


Compound 

A  GJ 
(kcal/mol) 

E'(SCF)b 

(hartrees) 

E(SCF)C 

(hartrees) 

AEg'd 

(kcal/mol) 

A  °Gsol 
(kcal/mol) 

DMA 

-5.37 

-135.020237 

-135.020398 

0.10 

-5.27 

TMA 

-3.53 

-174.281905 

-174.281973 

0.04 

-3.49 

MP 

-4.59 

-251.623022 

-251.623115 

0.11 

-4.48 

“Value  calculated  for  the  confonner’s  PCM-mediated  PBEPBE/6-31+G(d,p)  optimized  geometry. 
bPBEPBE/6-3 1 1+G(2d,2p)  single  point  energy  calculated  for  the  PCM-mediated  PBEPBE/6-3 1+G(d,p)  optimized 
geometry. 

cPBEPBE/6-31 1+G(2d,2p)  single  point  energy  calculated  for  the  gas-phase  PBEPBE/6-3 1+G(d,p)  optimized 
geometry. 
dE'(SCF)  -  E(SCF). 


'McQuaid,  M.  J.  The  Structure  of  Secondary  2-Azidoethanamines:  A  Hypergolic  Fuel  vs.  a  Nonhypergolic  Fuel;  ARL-TR- 
3176;  U.S.  Anny  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  2004. 

2 

“McQuaid,  M.  J.;  McNesby,  K.  L.;  Rice,  B.  M.;  Chabalowski,  C.  F.  Density  Functional  Theory  Characterization  of  the 
Structure  and  Gas-Phase,  Mid-Infrared  Absorption  Spectrum  of  2-Azido-N,N-dimethylethanamine  (DMAZ).  J.  Mol.  Struct. 
(TFIEOCFIEM)  2002,  587,  199-218. 
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Table  A-2.  DFT/PCM  results  for  DMAZ:  gas  phase. 


Label 

E'(SCF)a 

(hartrees) 

E(SCF)b 

(hartrees) 

AE'C 

(kcal/mol) 

ZPVE  +  A  Gcorr 
(hartrees) 

G, 

(hartrees) 

A  Gg 

(kcal/mol) 

Rel.  Pop. 

A 

-377.026882 

-377.027258 

0.24 

0.113168 

-376.914090 

0.5 

0.10 

B 

-377.028564 

-377.028763 

0.12 

0.113810 

-376.914953 

0.0 

0.25 

C 

-377.026792 

-377.026963 

0.11 

0.113516 

-376.913447 

0.9 

0.05 

D 

-377.026620 

-377.026849 

0.14 

0.112935 

-376.913914 

0.7 

0.08 

E 

-377.027261 

-377.027421 

0.10 

0.112617 

-376.914804 

0.1 

0.22 

F 

-377.027225 

-377.027356 

0.08 

0.113292 

-376.914064 

0.6 

0.10 

G 

-377.027066 

-377.027240 

0.11 

0.113032 

-376.914208 

0.5 

0.12 

H 

-377.025334 

-377.025546 

0.13 

0.113305 

-376.912241 

1.7 

0.01 

1 

-377.025188 

-377.025383 

0.12 

0.112667 

-376.912716 

1.4 

0.02 

J 

-377.024796 

-377.025100 

0.19 

0.112594 

-376.912506 

1.5 

0.02 

K 

-377.024801 

-377.025028 

0.14 

0.113205 

-376.911823 

2.0 

0.01 

L 

-377.025245 

-377.025455 

0.13 

0.113507 

-376.911948 

1.9 

0.01 

M 

-377.026490 

-377.028721 

1.40 

0.113993 

-376.914953d 

— 

— 

“PBEPBE/6-3 1 1+G  (2d,2p)  single  point  energy  calculated  for  the  confonner’s  PCM-mediated  PBEPBE/6-3 1+G(d,p)  optimized 
geometry. 

bPBEPBE/6-3 1 1+G(2d,2p)  single  point  energy  calculated  for  the  gas-phase  PBEPBE/6-3 1+G(d,p)  optimized  geometry. 

i i  t  sci  i  - 1 1  st  f  i. 

dConfonner  M  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  confonner  M 
ended  as  confonner  B. 


Table  A-3.  DFT/PCM  results  for  DMAZ:  aqueous  phase. 


Label 

A  Gj 
(kcal/mol) 

Gaq 

(kcal/mol) 

A  Gaq 

(kcal/mol) 

Rel.  Pop. 

A 

-7.83 

-7.05 

0.7 

0.12 

B 

-6.09 

-5.97 

1.7 

0.02 

C 

-6.86 

-5.81 

1.9 

0.02 

D 

-7.03 

-6.23 

1.5 

0.03 

E 

-7.56 

-7.37 

0.3 

0.21 

F 

-7.38 

-6.74 

1.0 

0.07 

G 

-7.44 

-6.86 

0.8 

0.09 

H 

-7.45 

-5.61 

2.1 

0.01 

I 

-7.70 

-6.17 

1.5 

0.03 

J 

-7.13 

-5.40 

2.3 

0.01 

K 

-6.67 

-4.56 

3.1 

0.00 

L 

-6.97 

-4.95 

2.8 

0.00 

M 

-9.11 

-7.71 

0.0 

0.38 

“Value  calculated  for  the  confonner’s  PCM-mediated  PBEPBE/6- 
31+G(d,p)  optimized  geometry. 
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Table  A-4.  DFT/PCM  results  for  MMAZ:  gas  phase. 


Label 

E'(SCF)a 

(hartrees) 

E(SCF)b 

(hartrees) 

AE'C 

(kcal/mol) 

ZPVE+AG™" 

(hartrees) 

G, 

(hartrees) 

AG, 

(kcal/mol) 

Rel.  Pop. 

A1 

-337.766635 

-337.766999 

0.23 

0.087913 

-337.679086 

1.9 

0.02 

B1 

-337.768361 

-337.768632 

0.17 

0.088562 

-337.680070 

1.2 

0.04 

Cl 

-337.769138 

-337.769566 

0.27 

0.088267 

-337.681299 

0.5 

0.16 

D1 

-337.76988 

-337.770222 

0.21 

0.088183 

-337.682039 

0.0 

0.35 

El 

-337.767703 

-337.767918 

0.13 

0.087767 

-337.680151 

1.2 

0.05 

FI 

-337.767613 

-337.767813 

0.13 

0.088031 

-337.679782 

1.4 

0.03 

G1 

-337.767475 

-337.767697 

0.14 

0.088021 

-337.679676 

1.5 

0.03 

HI 

-337.766140 

-337.766379 

0.15 

0.088038 

-337.678341 

2.3 

0.01 

11 

-337.765955 

-337.766183 

0.14 

0.087350 

-337.678833 

2.0 

0.01 

J1 

-337.768213 

-337.768454 

0.15 

0.087659 

-337.680795 

0.8 

0.09 

K1 

-337.767536 

-337.767749 

0.13 

0.088607 

-337.679142 

1.8 

0.02 

LI 

-337.767361 

-337.767609 

0.16 

0.087901 

-337.679708 

1.5 

0.03 

Ml 

-337.766406 

-337.768642 

1.40 

0.088598 

-337.680044d 

— 

— 

A2 

-337.766142 

-337.766437 

0.19 

0.088331 

-337.678106 

2.5 

0.01 

B2 

-337.768022 

-337.768331 

0.19 

0.088391 

-337.679940 

1.3 

0.04 

C2 

-337.766403 

-337.766632 

0.14 

0.088425 

-337.678207 

2.4 

0.01 

D2 

-337.766414 

-337.766694 

0.18 

0.087738 

-337.678956 

1.9 

0.01 

E2 

-337.767048 

-337.767278 

0.14 

0.087450 

-337.679828 

1.4 

0.03 

F2 

-337.766951 

-337.767185 

0.15 

0.088001 

-337.679184 

1.8 

0.02 

G2 

-337.766828 

-337.767048 

0.14 

0.087965 

-337.679083 

1.9 

0.02 

H2 

-337.766048 

-337.766290 

0.15 

0.087966 

-337.678324 

2.3 

0.01 

J2 

-337.764928 

-337.765226 

0.19 

0.087378 

-337.677848 

2.6 

0.00 

K2 

-337.764766 

-337.765033 

0.17 

0.087842 

-337.677191 

3.0 

0.00 

L2 

-337.765305 

-337.765580 

0.17 

0.088199 

-337.677381 

2.9 

0.00 

M2 

-337.765940 

-337.768336 

1.50 

0.088470 

-337. 679866s 

— 

— 

X2 

-337.766289 

-337.766525 

0.15 

0.087794 

-337.678731 

2.1 

0.01 

“PBEPBE/6-31 1+G(2d,2p)  single  point  energy  calculated  for  the  confonner’s  PCM-mediated  PBEPBE/6-31+G(d,p)  optimized 
geometry. 

bPBEPBE/6-3 1 1+G(2d,2p)  single  point  energy  calculated  for  the  gas-phase  PBEPBE/6-3 1+G(d,p)  optimized  geometry. 
cE'(SCF)  -  E(SCF). 

dConformer  Ml  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  conformer 
M 1  ended  as  confonner  B 1 . 

eConformer  M2  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  confonner 
M2  ended  as  conformer  B2. 
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Table  A-5.  DFT/PCM  results  for  MMAZ:  aqueous  phase. 


Label 

A  GJ 
(kcal/mol) 

Gaq 

(kcal/mol) 

A  Gaq 

(kcal/mol) 

Rel.  Pop. 

A1 

-10.38 

-8.3 

1.8 

0.02 

B1 

-9.17 

-7.8 

2.4 

0.01 

Cl 

-9.23 

-8.5 

1.6 

0.03 

D1 

-8.86 

-8.6 

1.5 

0.04 

El 

-9.79 

-8.5 

1.7 

0.03 

FI 

-9.77 

-8.2 

1.9 

0.02 

G1 

-9.67 

-8.0 

2.1 

0.02 

HI 

-10.05 

-7.6 

2.5 

0.01 

11 

-10.13 

-8.0 

2.2 

0.01 

J1 

-9.13 

-8.2 

1.9 

0.02 

K1 

-9.00 

-7.0 

3.1 

0.00 

LI 

-9.44 

-7.8 

2.3 

0.01 

Ml 

-11.53 

-10.1 

0.0 

0.53 

A2 

-10.04 

-7.4 

2.7 

0.01 

B2 

-8.79 

-7.3 

2.8 

0.00 

C2 

-8.98 

-6.4 

3.7 

0.00 

D2 

-8.94 

-6.8 

3.3 

0.00 

E2 

-9.52 

-8.0 

2.1 

0.01 

F2 

-9.39 

-7.5 

2.7 

0.01 

G2 

-9.48 

-7.5 

2.6 

0.01 

H2 

-9.91 

-7.4 

2.7 

0.01 

J2 

-9.53 

-6.7 

3.4 

0.00 

K2 

-9.20 

-6.0 

4.1 

0.00 

L2 

-9.34 

-6.2 

3.9 

0.00 

M2 

-11.01 

-9.5 

0.6 

0.19 

X2 

-8.47 

-6.2 

3.9 

0.00 

aValue  calculated  for  the  conformer’s  PCM-mediated  PBEPBE 
/6-31+G(d,p)  optimized  geometry. 
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Table  A-6.  DFT/PCM  results  for  PyrrAZ:  gas  phase. 


Label 

E'(SCF)a 

(hartrees) 

E(SCF)b 

(hartrees) 

AE'C 

(kcal/mol) 

ZPVE  +  A  Gcorr 
(hartrees) 

(hartrees) 

AG, 

(kcal/mol) 

Rel.  Pop. 

A 

-454.368672 

-454.369013 

0.21 

0.146784 

-454.222229 

0.67 

0.08 

B 

-454.370402 

-454.370566 

0.10 

0.147266 

-454.223300 

0.00 

0.24 

C 

-454.368870 

-454.369020 

0.09 

0.147070 

-454.221950 

0.85 

0.06 

D 

-454.368702 

-454.368884 

0.11 

0.146668 

-454.222216 

0.68 

0.08 

E 

-454.369200 

-454.369322 

0.08 

0.146252 

-454.223070 

0.14 

0.19 

F 

-454.369121 

-454.369240 

0.07 

0.146809 

-454.222431 

0.55 

0.10 

G 

-454.368996 

-454.369121 

0.08 

0.146742 

-454.222379 

0.58 

0.09 

H 

-454.365945 

-454.366099 

0.10 

0.146643 

-454.219456 

2.41 

0.00 

I 

-454.365807 

-454.365960 

0.10 

0.145899 

-454.220061 

2.03 

0.01 

J 

-454.365496 

-454.365709 

0.13 

0.145378 

-454.220331 

1.86 

0.01 

K 

-454.365329 

-454.365466 

0.09 

0.146219 

-454.219247 

2.54 

0.00 

L 

-454.365919 

-454.366088 

0.11 

0.146644 

-454.219444 

2.42 

0.00 

M 

-454.368351 

-454.368571 

0.14 

0.145866 

-454.222705 

0.37 

0.13 

A1 

-454.363064 

-454.363474 

0.26 

0.146238 

-454.217236 

3.80 

0.00 

B1 

-454.365580 

-454.365773 

0.12 

0.147418 

-454.218355 

3.10 

0.00 

Cl 

-454.363484 

-454.363673 

0.12 

0.147328 

-454.216345 

4.36 

0.00 

D1 

-454.363938 

-454.364137 

0.12 

0.148302 

-454.215835 

4.68 

0.00 

El 

-454.364485 

-454.364662 

0.11 

0.146523 

-454.218139 

3.24 

0.00 

FI 

-454.364457 

-454.364580 

0.08 

0.146978 

-454.217602 

3.58 

0.00 

G1 

-454.364391 

-454.364533 

0.09 

0.147017 

-454.217516 

3.63 

0.00 

HI 

-454.361856 

-454.362043 

0.12 

0.147184 

-454.214859 

5.30 

0.00 

11 

-454.361581 

-454.361746 

0.10 

0.146372 

-454.215374 

4.97 

0.00 

J1 

-454.360743 

-454.361401 

0.41 

0.146609 

-454.214792 

5.34 

0.00 

LI 

-454.361066 

-454.361260 

0.12 

0.147300 

-454.213960 

5.86 

0.00 

Ml 

-454.363159 

-454.365789 

1.65 

0.147305 

-454.218484d 

— 

— 

“PBEPBE/6-3 1 1+G  (2d,2p)  single  point  energy  calculated  for  the  conformer’s  PCM-mediated  PBEPBE/6-3 1+G(d,p)  optimized 
geometry. 

bPBEPBE/6-31 1+G(2d,2p)  single  point  energy  calculated  for  the  gas-phase  PBEPBE/6-31+G(d,p)  optimized  geometry. 
cE'(SCF)  -  E(SCF). 

dConformer  Ml  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  confonner  Ml 
ended  as  conformer  B 1 . 
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Table  A-7.  DFT/PCM  results  for  PyrrAZ:  aqueous  phase. 


Label 

A  GJ 
(kcal/mol) 

Gaq 

(kcal/mol) 

A  Gaq 

(kcal/mol) 

Rel.  Pop. 

A 

-8.76 

-7.9 

1.4 

0.06 

B 

-6.85 

-6.7 

2.6 

0.01 

C 

-7.78 

-6.8 

2.5 

0.01 

D 

-7.93 

-7.1 

2.2 

0.02 

E 

-8.55 

-8.3 

1.0 

0.12 

F 

-8.40 

-7.8 

1.5 

0.05 

G 

-8.43 

-7.8 

1.5 

0.05 

H 

-8.49 

-6.0 

3.3 

0.00 

1 

-8.65 

-6.5 

2.8 

0.01 

J 

-8.09 

-6.1 

3.2 

0.00 

K 

-7.38 

-4.8 

4.6 

0.00 

L 

-7.81 

-5.3 

4.0 

0.00 

M 

-9.83 

-9.3 

0.0 

0.66 

A1 

-10.06 

-6.0 

3.3 

0.00 

B1 

-7.49 

-4.3 

5.1 

0.00 

Cl 

-8.12 

-3.6 

5.7 

0.00 

D1 

-8.05 

-3.2 

6.1 

0.00 

El 

-8.92 

-5.6 

3.7 

0.00 

FI 

-8.64 

-5.0 

4.3 

0.00 

G1 

-8.82 

-5.1 

4.2 

0.00 

HI 

-8.92 

-3.5 

5.8 

0.00 

11 

-9.15 

-4.1 

5.2 

0.00 

J1 

-7.70 

-1.9 

7.4 

0.00 

LI 

-7.45 

-1.5 

7.9 

0.00 

Ml 

-10.66 

-6.0 

3.3 

0.00 

“Value  calculated  for  the  conformer’s  PCM-mediated  PBEPBE 
/6-31+G(d,p)  optimized  geometry. 
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Table  A-8.  DFT/PCM  results  for  MorphAZ:  gas  phase. 


Label 

E'(SCF)a 

(hartrees) 

E(SCF)b 

(hartrees) 

AE'C 

(kcal/mol) 

ZPVE  +  A  Gcorr 
(hartrees) 

G* 

(hartrees) 

A  Gg 

(kcal/mol) 

Rel.  Pop. 

A 

-529.532886 

-529.533193 

0.2 

0.151273 

-529.381920 

0.6 

0.11 

B 

-529.534570 

-529.534774 

0.1 

0.152011 

-529.382884 

0.1 

0.26 

C 

-529.532752 

-529.532956 

0.1 

0.151793 

-529.381163 

1.1 

0.05 

D 

-529.532625 

-529.532882 

0.2 

0.151145 

-529.381737 

0.7 

0.09 

E 

-529.533036 

-529.533223 

0.1 

0.150741 

-529.382482 

0.3 

0.19 

F 

-529.532920 

-529.533107 

0.1 

0.151044 

-529.382063 

0.5 

0.12 

G 

-529.532833 

-529.533016 

0.1 

0.151144 

-529.381872 

0.6 

0.10 

H 

-529.531035 

-529.531242 

0.1 

0.151338 

-529.379904 

1.9 

0.01 

I 

-529.530935 

-529.531140 

0.1 

0.150700 

-529.380440 

1.5 

0.02 

J 

-529.530547 

-529.530833 

0.2 

0.150674 

-529.380159 

1.7 

0.02 

K 

-529.530147 

-529.530378 

0.1 

0.150379 

-529.379999 

1.8 

0.01 

L 

-529.530905 

-529.531105 

0.1 

0.151577 

-529.379528 

2.1 

0.01 

M 

-529.532530 

-529.534773 

1.4 

0.151889 

-529.382884d 

— 

— 

A1 

-529.526172 

-529.526592 

0.3 

0.151008 

-529.375584 

4.6 

0.00 

B1 

-529.528559 

-529.528860 

0.2 

0.151958 

-529.376902 

3.8 

0.00 

Cl 

-529.526423 

-529.526653 

0.1 

0.150994 

-529.375659 

4.5 

0.00 

D1 

-529.526725 

-529.526976 

0.2 

0.150993 

-529.375983 

4.3 

0.00 

El 

-529.527316 

-529.527515 

0.1 

0.150968 

-529.376547 

4.0 

0.00 

FI 

-529.527224 

-529.527407 

0.1 

0.151385 

-529.376022 

4.3 

0.00 

G1 

-529.527139 

-529.527329 

0.1 

0.151341 

-529.375988 

4.3 

0.00 

HI 

-529.525027 

-529.525230 

0.1 

0.151455 

-529.373775 

5.7 

0.00 

11 

-529.524820 

-529.525038 

0.1 

0.150653 

-529.374385 

5.3 

0.00 

J1 

-529.524570 

-529.524804 

0.1 

0.150977 

-529.373827 

5.7 

0.00 

LI 

-529.524401 

-529.524609 

0.1 

0.150790 

-529.373819 

5.7 

0.00 

Ml 

-529.526256 

-529.528811 

1.6 

0.151591 

-529.377220" 

— 

— 

“PBEPBE/6-3 1 1+G  (2d,2p)  single  point  energy  calculated  for  the  conformer’s  PCM-mediated  PBEPBE/6-3 1+G(d,p)  optimized 
geometry. 

bPBEPBE/6-311+G(2d,2p)  single  point  energy  calculated  for  the  gas-phase  PBEPBE/6-3 1+G(d,p)  optimized  geometry. 

'l  '(SC  I  t  -  E(SCF). 

dConformer  M  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  conformer  M 
ended  as  conformer  B. 

eConfonner  Ml  is  an  aqueous-phase  structure  that  has  no  gas-phase  analog.  Gas-phase  optimizations  started  from  conformer  Ml 
ended  as  conformer  B 1 . 
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Table  A-9.  DFT/PCM  results  for  MorphAZ:  aqueous  phase. 


Label 

A  GJ 
(kcal/mol) 

Gaq 

(kcal/mol) 

A  Gaq 

(kcal/mol) 

Rel.  Pop. 

A 

-13.09 

-12.3 

0.9 

0.07 

B 

-11.68 

-11.5 

1.7 

0.02 

C 

-12.24 

-11.0 

2.1 

0.01 

D 

-12.44 

-11.6 

1.6 

0.02 

E 

-13.42 

-13.1 

0.1 

0.25 

F 

-13.33 

-12.7 

0.5 

0.14 

G 

-13.31 

-12.6 

0.6 

0.11 

H 

-13.43 

-11.4 

1.7 

0.02 

I 

-13.5 

-11.8 

1.3 

0.03 

J 

-12.78 

-10.9 

2.3 

0.01 

K 

-12.28 

-10.3 

2.8 

0.00 

L 

-12.81 

-10.6 

2.6 

0.00 

M 

-14.58 

-13.2 

0.0 

0.31 

A1 

-14.46 

-9.6 

3.6 

0.00 

B1 

-12.29 

-8.3 

4.8 

0.00 

Cl 

-12.63 

-8.0 

5.2 

0.00 

D1 

-12.63 

-8.1 

5.0 

0.00 

El 

-13.73 

-9.6 

3.5 

0.00 

FI 

-13.52 

-9.1 

4.1 

0.00 

G1 

-13.68 

-9.2 

3.9 

0.00 

HI 

-13.44 

-7.6 

5.6 

0.00 

11 

-13.25 

-7.8 

5.4 

0.00 

J1 

-11.80 

-6.0 

7.2 

0.00 

LI 

-11.84 

-6.0 

7.2 

0.00 

Ml 

-15.32 

-10.2 

3.0 

0.00 

aValue  calculated  for  the  conformer’s  PCM-mediated  PBEPBE 
/6-31+G(d,p)  optimized  geometry. 
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AMSRD  ARL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  PE 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 
1  DIR  USARL 

AMSRD  ARL  Cl  OK  TP  (BLDG  4600) 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  ARMY  RSRCH  OFC 

D  MANN 
POBOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
22709-2211 

1  ARMY  RSRCH  OFC 
CHEMICAL  SCI  DIV 
R  ANTHENIEN 
POBOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
22709-2211 

7  US  ARMY  AVN  &  MIS  CMND 
AMSRD  AMR  PS  PT 
G DRAKE 
L  FELTON 
J  LILLY 
N  MATHIS 
R  MICHAELS 
J  NEIDERT 
D  THOMPSON 
BLDG  7120 
REDSTONE  RD 

REDSTONE  ARSENAL  AL  35898 

2  HUGHES  ASSOC  INC 
W  RUPPERT 

N  LIEB 

3610  COMMERCE  DR  STE817 
BALTIMORE  MD  21227 

1  3M  COTTAGE  GROVE  CTR 
70-2SE 

A  MANZARA 

COTTAGE  GROVE  MN  55016 

1  USEPA  HDQTRS 
J  NABHOLZ 
ARIEL  RIOS  BLDG 
1200  PENNSYLVANIA  AVE  NW 
MC  7403M 

WASHINGTON  DC  20460 


ABERDEEN  PROVING  GROUND 

1  COMMANDER 

US  ARMY  RDECOM 
AMSRD  MSF 
E  HANGELAND 
5183  BLACKHAWK  RD 
E-5101  RM  182 
APG-EAMD  21010-5424 


1  DIRECTOR 
USACHPPM 

DIRECTOR  OF  TOXICOLOGY 
MCHB  TS  THE 
M  JOHNSON 
APG-EAMD  21010-5403 

19  DIRUSARL 

AMSRD  ARL  WM  B 
M  ZOLTOSKI 
AMSRD  ARL  WM  BD 
SBUNTE 
E  BYRD 
C  CHEN 
B  FORCH 
M  HURLEY 
A  KOTLAR 
PKASTE 
K  MCNESBY 
J  MORRIS 
M  NUSCA 

R  PESCE-RODRIGUEZ 
S  PIRAINO 
BRICE 
R  SAUSA 

AMSRD  ARL  WM  M 
J  BEATTY 
S  MCKNIGHT 
AMSRD  ARL  WM  MA 
M  VANLANDINGHAM 
AMSRD  ARL  WM  T 
P  BAKER 
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